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Abstract 

We consider a left invariant Riemannian metric on SO 3 with two equal eigenval¬ 
ues. We find the cut locus and the equation for the cut time. We find the diameter 
of such metric and describe the set of all most distant points from the identity. Also 
we prove that the cut locus and the cut time converge to the cut locus and the cut 
time in the sub-Riemannian problem on SO 3 as one of the metric eigenvalues tends 
to infinity. 
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1 Introduction 

Consider a left invariant Riemannian metric on SO3. Restriction of the metric to the 
tangent space at the identity is diagonalizable with respect to the Killing form. Let R, 
J2, I3 be its eigenvalnes. In this paper we consider only the Lagrange case (R = R). 

Parametrization of left invariant Riemannian geodesics on SO3 is a classical L. Enler’s 
resnlt [T]. L. Bates and F. Fasso [ 2 ] fonnd an eqnation for the conjngate time (in the 
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Lagrange case) and the conjngate locns depending on the ratio of eigenvalnes. Thns 
local optimality of geodesics was stndied. Bnt global optimality of geodesics was not 
investigated yet. 

For the same problem on SU2, T. Sakai [ 3 ] proved that in the case Ji > J3 the cnt locns 
is a two dimensional disk. In the case Ji < Js there is a conjectnre stated in M. Berger’s 
book |3] that the cnt locns is a segment. 

Notice that if Ii = I2 = I3 (the Enler case) a geodesic on SO3 is a set of rotations 
aronnd a hxed line. Snch a geodesic is optimal when the rotation angle is in the segment 
[0,7r]. So in the Enler case the cnt locns is the set of all axial symmetries, and the cnt 
time eqnals 7r\/Ji. 

In this work we stndy the global optimality of geodesics on SO3 in the Lagrange case, 
describe the Maxwell strata and the cnt locns, give an eqnation for the cnt time. We nse 
the method introdnced by the second co-anthor for the generalized Dido problem [ 5 ] and 
the problem of Enler elasticae [6]. First we describe a symmetry gronp of an exponential 
map. Secondly we hnd its hxed points (Maxwell strata). And then we prove that the 
exponential map is a diffeomorphism of some open set to the complement of the Maxwell 
strata in SO3. 

If 2/1 ^ J3 then the cut locus is a projective plane (the set of all axial symmetries). If 
2/1 < I3 then the cut locus is the union of that projective plane and a segment which is 
the set of certain rotations around the axis corresponding to the eigenvalue I3. Also we 
hnd the diameter of SO3. 

Besides we prove Berger’s conjecture about the cut locus on SU2 in the case A <13. 

Moreover, we show that the parametrization of geodesics, the conjugate time, the 
conjugate locus, the cut time and the cut locus converge (as I3 tends to inhnity) to the 
same objects in a sub-Riemannian problem on SO3. The sub-Riemannian problem on 
SO3 was hrst studied by V. N. Berestovskiy and 1 . A. Zubareva [H] (they described 
singularities of spheres), and then by U. Boscain and F. Rossi [ 7 ] (they got the equation 
for the cut time and described the cut locus). 

This paper has the following structure. The problem is stated in Section [ 2 l In Section [ 3 ] 
we present equations for geodesics in terms of quaternions (i.e., on a double covering of 
SO3). In Section m we describe symmetries of the exponential map. Next, in Section [ 5 ] 
we hnd the Maxwell strata corresponding to the symmetries of the exponential map 
fSubsection l 5 .ll) . then we study a relative location of the Maxwell strata fSubsection 15 . 31 ) . 
and we hnd the hrst Maxwell time. It turns out that the hrst Maxwell time corresponding 
to the symmetries is in fact the cut time. The main theorem about the cut time and the 
cut locus is stated in Section [3 In Section [6] some necessary results about the conjugate 
time are stated. In Section [8] we compute the diameter of SO3 and hnd the set of all most 
distant points from the identity. In Sections O and HH] we discuss some connections of our 
work to the same problem on SU2 and to the sub-Riemannian problem on SO3. 
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2 Left invariant Riemannian problem on SO3 


Any left invariant Riemannian metric on SO3 is completely determined by a positive 
definite quadratic form J on the tangent space TidSOa = 503. Let ei, 62, 63 be an 
orthonormal basis (with respect to the Killing form) such that J is diagonal in this basis. 
Let Ji, J2, Iz be the corresponding eigenvalues of J. 

The problem of hnding length minimizers is equivalent to the following optimal control 
problem: 

Q = QVt, Vt = uiei + ^262 + M3e3 e SO3, 

Q G SO3, (mi,M2,'U3) G 
( 5(0) = id, Q{ti) = Qi, 

{Iiu\ + /2M2 + dt —)■ min . 

Minimization of this functional is equivalent to minimization of arc length due to the 
Cauchy-Schwartz inequality. 

If there is a triangle with the edges of length Ji, I2, 13 then this problem has a mechan¬ 
ical interpretation: rotation of a rigid body around a hxed point. The numbers Ii, 12,13 
are the inertia moments of this rigid body [T] . 

Let us use the Pontryagin maximum principle mu- The Hamiltonian of the maximum 
principle is 

Huip) = uiPi + U2P2 + U3P3 hul /3M3), 

where u ^ 0 , p = piSi +P2S2 G SO3, {et} is a basis of SO3 dual to {cj} with respect 

to the Killing form. If u{t) is an optimal control then almost everywhere (pit)) = 
ma.Xu H^{p{t)) . If z/ = 0 then p = 0 , which contradicts the nontriviality condition of the 
maximum principle (the pair (z^,p) must be nonzero). Thus u < 0 . The pair (z^,p) is 
dehned up to a positive factor, let u = —1. Then 



Pi ■ 1 o Q 

Ui = —, t = I, 2, 3. 


(Here and below we identify SO3 and SO3 via the Killing form.) Then the maximized 
Hamiltonian is 

n 2 \h h hj 

The Hamiltonian system with the Hamiltonian H reads 


Q = Qn, 
p = (ad* fl)p, 

where {Q,p) G SO3 x SO3 (trivialization of the cotangent bundle T*S03 via the left action 
of SO3), and H = ^ei -f ^62 + 7^63 G SO3. Notice that L. Euler got the same equations 
from the conservation laws. 
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3 Solution of the Hamiltonian system 
in the Lagrange case 


Below we consider only the Lagrange case Ii = I2 (the case Ii = I2 = I3 is called the 
Euler case). In the Lagrange case the Hamiltonian system ([T]) is integrated in elementary 
functions [2]: 

Q{t) ( 2 ) 

P{t) = P, ( 3 ) 

where p(0) = piCi +^>262 +^363 G SO3, Q(0) = id, denotes the rotation of by an 
angle if aronnd a vector v G SO3 (the direction of the rotation is snch that for any vector 
w ^ span{n} the frame {w, Ry^^w,v) is positively oriented). The parameter 
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dehnes oblateness for the rigid body. We identify elements of SO3 with orthogonal trans¬ 
formations of SO3 by the co-adjoint representation. 

It is convenient to make further computations in terms of quaternions. Consider the 
double covering 

n : {g G e I |g| = 1} ~ ^ SO3. 

Any quaternion q of the unit norm can be written in the form 


q = cos 



-|- sin 



{aR -f- 02j ask), 


where oi, 02, 03 G M, -l- a| = 1. By dehnition n(±g) = Ra^ip is the rotation by the 

angle (p around the vector a = OiCi -|- 0262 -l- 0363. 

Let ±(go('r) + gi('r)i -l- g2('r)j + Q3{T)k) G El be the lift to of the solution ([ 2 ]) of the 
horizontal part of the Hamiltonian system ([T]), where we change the time scale t = ^\p\- 
Multiplying the two quaternions corresponding to the two factors in the solution Q{t), 
we obtain: 

{ go(r) = cos(r) cosirpps) - ps sin(r) sinirpps), 

( ll[r] ) ^ sin(r)i?e3,-r^p-3 ( ) ’ ( 4 ) 

qsir) = cos(r)sin(rr7P3)-Fp3sin(r)cos(r?7P3), 

where p = p, and we denote the restriction of Res,a to the plane span{ei, 62} by the same 
symbol. 

Solution ([ 3 ]) of the vertical part of the Hamiltonian system ([T]) reads 


p(r) = i?e3 -2r,?p3P(0), IpI = const. 


( 5 ) 
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4 Symmetries of the exponential map 

Consider geodesics which start at the identity and have the nnit velocity. In this section 
we dehne some symmetries on the set of such geodesics. 

Let u G TidSOs be a tangent vector to such a geodesic at the identity. Then 

J{u) = Iiu\ + J2M2 + = 1 . 

By the maximum principle we have Ui = ^, i = 1 , 2 , 3 . Then 

m = d + d + d = 2//(p). 

^2 ^3 

When we identify 503 and SO3, the unit sphere in TidSOa corresponds to the level surface 
of the Hamiltonian (an ellipsoid) 

C = {peso;\ H{p) = 1 / 2 }. 

Recall that the exponential map is a map 

Exp : C X M+ —)■ SO3, Exp{p,t) = tt o e*^(id,p), 

where p E C, t E M+, and is the flow of the Hamiltonian vector field H, while 
TT : T*S03 —)■ SO3 is the projection of the cotangent bundle to its base. 

Definition 1. A pair of diffeomorphisms s : C x M+ -E C x M+ and 's : SO3 —)■ SO3 is 
called a symmetry of the exponential map if 

Exp o s = s'o Exp. 


Denote the vertical part of the Hamiltonian vector held 

VP 3 


Hyert (p) 


h^pl+pl 


(-^261+^162). 


We will consider only the symmetry group S of the exponential map that corresponds to 
isometries of C that preserve Hyert or change the sign of the vector held Hy^rt- 

The generators of the group S are the rotations around 63, the rehection ai in the 
plane span{ei,e3} and the rehection <72 in the plane span{ei,e2}. It is easy to see that 
there is an isomorphism S ^ O2 x Z2. 

As we identify SO3 with the set of orthogonal transformations of SO3 by the co-adjoint 
representation, we can say that S' is a subgroup of O3. 
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Proposition 1 . The group S is embedded into the group of symmetries of the exponential 
map. An element a G S' maps to the pair of diffeomorphisms 

s„ : C X M+ —)■ C X M+ and % : SO3 —)■ SO3, 

which are defined as follows: 


SaiP^t) = 


if (J^Hyert = -HyerU 


e SO3. 


Proof. It is clear that the defined map a 1—)■ {sa,'Sa) is an injection of the group S to the 
direct product of the diffeomorphism groups of C x R+ and SO3. So we need to show 
that for any generator a G S the pair of diffeomorphisms is a symmetry of the 

exponential map. 

Notice that the action = Ra{v),^ means that a acts on the imaginary part of 

the quaternion that corresponds to G SO3. 

Let Exp(p,t) = n(go(t) + qi{t)i + q2{t)j + q‘i{t)k), Exp o s„{j)R) = n(go(t) + qi{t)i + 
q2{t)i + q^{t)k). 

Let us use (jl]) to compute the action of a symmetry on the image of the exponential 
map. 

( 1 ) If cr = i?e3,a then %{t) = qo{T), %{t) = g3(r), and 


( Ur) 

V 52(r) 


Sill(^T R, 


63,0 


Pi 

P 2 


Re^,a sin^T ) Rg^^—rrip^ 


Pi 

P2 


We have Exp o sr^^^^{p, t) = o Exp(p, t). 

( 2 ) If a is the reflection in the plane span{ei, 63} or in the plane span{ei, 62} then it 
is an isometry that changes the sign of the vertical part of the Hamiltonian vector held. 
We get Exp o s^ipR) = Exp((Te*^“"’"*(p), t) = {aRe3-2Tr,aip)3{p)R)- 

( 2 a) If a = (Jl is the rehection in the plane span{ei,e3} then Ur) = qo{r), Ur) = 
qsir). Using cri(p)3 = ps, we obtain 

^ = sin(r)i?e3,-T»?p30'l.Re3,-2r?7P3 ^ = Cn sin(r)i?e3,-r77P3 

where we denote the restriction of ai to the plane span{ei,e2} by the same symbol. We 
have Exp o (p, t) = U ° Exp(p, t). 



( Ur) 
V Ur) 
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(2b) Let (T = (72 be the reflection in the plane span{ei,e2}. Since go and are 
respectively even and odd functions of the variable ps we see that go(T) = go(T), %{t) = 
—<?3 ('?■)• We have 


Ur) 

Ur) 


sin(T ) -Res, —Tr)(72(p)3 ^ 2 -^ 63 ,— 2 rJ 7 P 3 


Pi 

P 2 


From (72 (p) 3 = 

f Ur) \ 
V Ur) J 

We have Exp 


= —p3, we obtain 
sin (t) Res,rr^psR 
O (P, t) = U ° 


e3,-2T?7p3 ^ 

Exp(pR). 


Pi 

P2 


sin(T )Re3,— T77P3 


Pi 

P2 


Ur) \ 
Ur) J ' 

□ 


5 Maxwell strata 

Recall that we assume that all geodesics have an arc length parametrization. 

Definition 2 . A point Q G SO3 is called a Maxwell point if there exist two different 
geodesics Qi,(32 : [ 0 ;^] SO3 such that they reach the point Q at the same time 
Q = QiitraU = Q 2 (tma.x)- This time is called a Maxwell time. 

It is well known that a geodesic is not optimal after a Maxwell point. 

Definition 3 . The first Maxwell set in the pre-image of the exponential map is the set 

M = {(pRmax) e C X M+ I V e C \ {p} : Exp(p,fmax) = Exp(p', f^ax), 

but Vt e (ORmax) Vpi e C \ {p} Exp(pR) R Exp(piR)}. 

The time fmax is called the first Maxwell time for p E C. 

Obviously, any point of Exp Ad is a Maxwell point. 

Definition 4 . Let A be a subset of the group S. The first Maxwell set corresponding to 
the set A in the pre-image of the exponential map is the set 

M{A) = {(p, fmax) e O X M+ I 3(7 e A, a R id : Exp(p, f^ax) = Exp o s„(p, f^ax), 

but Vt e ( 0 , fmax) Exp(p, t) R Exp o s^(p, t)}. 

The time fmax is called the first Maxwell time for p E C corresponding to the set of 
symmetries A. 

Generally speaking, this time can be greater than the first Maxwell time. 
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The aim of this section is a description of the hrst Maxwell sets in image and pre-image 
of the exponential map. First, for any a G S' we describe Ai{a). Secondly, we stndy the 
relative location of the sets A 4 (a), and then we find 

M{S) C \jM{a). 

<j^S 

Thirdly, we prove that the exponential map is a diffeomorphism of an open subset of 
C X R+ bounded by A^(S) onto SO3 \ (Exp A^(S) U {id}). This implies that A^(S) 
and Exp Ai{S) are the cut loci in the pre-image and the image of the exponential map, 
respectively. 

5.1 Maxwell strata corresponding to the symmetries 
of the exponential map 

Definition 5. Denote the smallest positive roots of the equations qo^r) = 0 and gsi^r) = 0 
by ro(p3) and nips), respectively. 

We consider tq and as functions of the variable p^. These functions depend on the 
parameter r]. If p^ = 0 , then the equation gsi^r) = 0 is identically satisfied, and the value 
r3(0) is not defined. So, we have 

D : [“1; 1] (0, +00], 

D : [—1,1] \ {0} (0, -|-cxo]. 

Proposition 2. The set 

At (cr) = Ato U At 12 U At3 

(tGS 

contains three components: 

Mo := {(p,t) eCxM+jt = 

IPI 

Mi2 := {(p,t) ec xR+ I P3 7^ ± 1 , t = 

At3 := {(p,t) e C X M+ I P3 7^ 0 , t = 

Let us consider elements of SO3 as transformations of the form where |n| = 1 
and cp G [0,7r]. Such transformations form a three dimensional projective space. We will 
realize this projective space as a three dimensional ball with antipodal identification of 
the boundary points (because = R-v,n)- 

To prove Proposition [ 2 ] we need the following lemma. 









Lemma 1. A geodesic Q{t) = Exp(p, t) reaches the circle 


I V -L 63, |u| = 1 } 


only for Ps = 0 . 

Proof. If Q{t) reaches this circle then there exists r such that 


9 o(t) = 0 , 

QsiT) = 0 . 


(6) 


If cos(r) cos^rrips) 7 ^ 0 then this system of equations is equivalent to the system of 
equations 

{ 1 — p 3 tan(r) tan{TT]p 3 ) = 0 , 

tan^rrjps) +p3tan(r) = 0. 

It follows that 1 +p|tan^(r) = 0, a contradiction. 

Hence, cos^rrip^) = 0 or cost = 0. Let us consider these two cases. 

In the first case from qsir) = 0 we have cost = 0, and from qo{T) = 0 we obtain 
P 3 sinr = 0. This yields that ps = 0 and cos^rrjps) = 1, a contradiction. 

In the second case cost = 0 we have t = |, and system of equations (jh]) reads 


TT _ 

-pssm{-r]p3) = 0 , 

TT _ 

P3COS{-pp3) = 0 . 


This implies that ps = 0 . □ 

Proof of Proposition 0 For any symmetry a G S' we will find its fixed points 

SO^ = {Q G SO3 I = Q}. 

Notice, that if G SO3 then v = a{v) or v = —cr{v) and ip = n. 

It is clear that Exp A 4 {a) C SO3. Let us study geodesics that are symmetric with 
respect to a and reach a fixed point of a at the same time fmax- We will find this time 
and describe j\ 4 (a). Let us denote Xmax = 

Elements of S can be of the following types: 

■^63,0;; ^ 1 -^ 63 ,a: ^2-^63,a: ^ 2 ^lPez,ct' 

Let US consider them consecutively. 

(la) Let a = Re^^a and a G ( 0 , tt). If Ry^^ G SO3 then Re3,a'v = v, i.e., v || 63. 
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Hence, for a quaternion that corresponds to G SO3, we have gi = ^2 = 0 and 
P1+P2 7^ 0 (otherwise Re3,aP = P and the symmetry preserves the corresponding geodesic). 
Thus, the hrst Maxwell time corresponding to this symmetry can be found from the 
equation sinxmax = 0 , consequently f ma x = It follows that M.[a) = M.i2- 

(lb) Let a = Re^,^- Then v || 63 or n T 63 and ip = n. The hrst situation was already 
described in (la), let us consider the second one. In this case by Lemma [ 1 ] it follows that 
p3 = 0 . 

So, we have the piece {{p,t) G M.q\ ps = 0 } of the component Mq. 

( 2 ) Let a = aiRe^^a, where Ui is the rehection in the plane span{ei, 63}. We can assume 
a = 0 and cr = cxi. In the general case the Maxwell set will be the result of rotation by 
the angle — ^ around 63 of the Maxwell set corresponding to the symmetry cxi. 

If Rv^ip = Ro-i{v),ifi, then v G span{ei, 63} or n T span{ei, 63} and ip = n. 

In the hrst case from the formulas of qi and q2 it follows that 

sin(Tinax).Re 3 ,Tinaxr?P 3 ( .pr ) ~ sin(Tniax)c’’l.R e 3 ,'rmax??P 3 ( ^ 

\ P 2 J \ P 2 

If siuTmax ^ 0 , then Re^,r^,^np^{p) G span{ei,e3} and i?e3,-2r™,w3^i(p) = P- Thus, 
'5(71 (p, ^max) = (P)^max), l-G., this Symmetry preserves the corresponding geodesic. 

If sinr„ia^ = 0 then t^ax = and p is such that ^ (p,tmax), i-e., 

Rei,T^axnPii.P) ^ span{ei,e3}. For diherent a we have the same component AI12 as in ( 1 ). 
In the second case we have the same situation as in (lb). 

( 3 ) Let a = o'2Re3,a where (J2 is the rehection in the plane span{ei, 62}. 

Then v = cr2Re3,a{v) or v = —o'2Re3,a{v) and ip = n. 

In the hrst case we have Res,a = id and v G span{ei, 62}. Hence, qs^Tmax) = 0 , and we 
have the component AI3 = {{p, | ^3 ^ 0} (when P3 = 0 the point p is preserved 

under the considered symmetry). 

In the second case, if n = ±63 then a can be arbitrary, otherwise a = tt. From ip = tt 
it follows that qo^Tmax) = 0 . We have the component M.q. 

( 4 ) Let cr = o'2(TiRe3,a = -Rei,77-^63,a he the axial symmetry in i?e3^_|(ei). We assume 
a = i?ei,7r (then we have to consider the results of rotations around 63 of the corresponding 
Maxwell set). 

There are two cases: n || ei or n T ei and ip = n. Thus, we have already found the 
components A 4 o and □ 

5.2 Continuity of the functions tq and ra 

To describe a relative location of the sets Ato, Aiu and Ms we need to compare 

and Recall that tq and Ts are respectively the minimal positive zeros of the 



10 





functions and ^3. The value \p\ depends only on p3, it follows that we need to compare 
^ and T^^ps), where G [— 1 , 1 ]. First of all we will stndy some properties of the 
fnnctions tq and T3. 

The graphs of the fnnctions tq and with different values of the parameter r] are 
shown in Fig. [TJ 


Figure 1 : The graphs of the functions tq and 




-| < ?7 < 0 


T] > 0 


Remark 1 . The fnnctions tq and T3 are even. Indeed, the fnnction qq is even, and the 
fnnction is odd relative to valne p^. Therefore the smallest positive zeros of these 
fnnctions do not depend on the sign of p^. 


Proposition 3. The functions tq and are continuous on the sets [—1,1] and [—1,1] \ 
{0}, respectively. 

Proof. It follows from the implicit fnnction theorem that it is enongh to prove that the 
functions go and q^ have no multiple zeros. 

Differentiation of the fnnctions go and go with respect to r gives 


dqo 

dr 


T 


— (1 + rjpl) sinr cos(rr7P3) — p^^rj + 1 ) cost sin(rr7P3), 


(7) 


1 . 

that 


dqs. . 

97 ^ = 

Assnme that go 


— (1 + ppI) sinrsin(r? 7 P 3 ) +^ 3(77 + 1 ) cost cos(rr^po). 
has a multiple zero, then there exists po G [—1,1] and 


( 8 ) 

r > 0 snch 


qo{^) = 0 , 


dqo 

dr 


(r) = 0. 


( 9 ) 
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Consider first the case cost cos (r?7P3) 7^ 0 . Dividing both eqnations by cos r cos (r77^3), 
we get 

{ 1 — p3 tanr tan(rr7P3) = 0, 

-(1 + r]pl) tanr - ^3(77 + 1) taii{Tr]p3) = 0. 

Note that tanr 7^ 0 . Expressing tan(r?7P3) from the hrst eqnation, we obtain 


tan^ r = 


p + l 

1 + ppl 


< 0 , 


we get a contradiction. 

Consider now the case cos r cos^rpp^) = 0 . 
If cos r = 0 then from ([ 9 |) we get 


( -p3sin(rr7P3) = 0, 

I -(1 + r]pl) cos{Tr]p3) = 0. 

It follows (since sine and cosine can not be eqnal to zero simnltaneonsly and 1 + pp^ > 0 ) 
that p3 = 0, hence cos^rpp^) = 1, a contradiction. 

If cosi^rpps) = 0 then from ([H]) we get 


-P3 sinr = 0, 

-psiV + 1 ) COST = 0 , 

thns p3 = 0 . Then cos{Tpp3) = 1 , a contradiction. 

2. Assnme that gs has a mnltiple zero, then there exists p3 E [— 1 , 1 ] \ {0} and r > 0 
snch that 

= 0 , 


dqs 

17 ' 


t) = 0 . 


( 10 ) 


Let cos T cosi^rpps) 7^ 0 . Dividing both eqnations by cost cos{Tpp3), we get 

tan(r77P3) + p 3 tan r = 0, 

-(1 + ppI) tanrtan(r?7P3) + psip + 1) = 0. 


This implies that 


9 P + 1 

tan r = —^^ < 0, 


l + hPs 

we get a contradiction. 

Now consider the case cos t cos^rpps) = 0 . 
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If COST = 0 then from fflU]) . we get 


( p3 cos(r?7P3) = 0, 

[ -(1 + r]pl) sin(r?7P3) = 0. 

This yields that p3 = 0 (since sine and cosine can not be equal to zero simultaneously and 

1 + pPI > 0 ). 

If cos{Tpp3) = 0 then from (fTOD . we obtain 

{ cos r = 0, 

-(1 + 77p2)sinr = 0, 

we get a contradiction. □ 

5.3 Relative location of the Maxwell strata 

To get a description of the relative location of the Maxwell strata we shall compare Tq^Ps), 
71 and for different values of p^ G [—1,1]. 

Proposition 4. For any P3 7 ^ 0 the inequality Tq^P^) < r^^ps) is satisfied. 

Proof. Notice that if p^ = ±1 then the statement is correct. Actually, in this case 

go(r) = cos(r(l + 77)), 

qsir) = ±sin(r(l + r])). 

Hence ro(±l) = < (ifii) = ^ 3 (± 1 )- 

Suppose that the inverse statement is satisfied, let p3 7^ 0 be such that Toi^ps) ^ 

From Proposition [ 3 ] we know that the functions Tq and T3 are continuous. It follows that 
there exists Ps 7^ 0 such that To(p3) = Tzifi^). This means that there exist p^ and r 
such that the geodesic, corresponding to the co-vector p, reaches the circle defined by the 
equations qo = qs = 0 . By Lemma [U this is possible only for p3 = 0 , a contradiction. □ 

The above proposition means that if p3 7^ 0 then geodesics reach Exp Aio earlier than 
Exp AT3. (If P3 = 0 then the corresponding geodesic always lies in Exp Af 3.) This means 
that Ais is not the first Maxwell set. 

Let us consider now the strata A 4 o and Ali2. 

Proposition 5. ( 1 ) If p ^ then for any p^ G [— 1 , 1 ] the inequality tq^Ps) ^ ti is 
satisfied. 

( 2 ) Ifp<-I then Toips) A vr /or ^ ^ jpg] ^ 1 and To^pfi) < tt for |p3| < 

See Fig. [21 
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Figure 2: The function Tq and vr. 



V <-I V =-I -i<? 7<0 77 = 0 , To = 1 ? 7>0 


Proof. Note that go(0) = 1- So, to prove (1) and the second part of (2) (i.e., that the 
function qoir) has a positive root less or eqnal than tt) it is enongh to find a number 
0(h)Ps) ^ ^ such that (?o(0(h)P3)) ^ 0- Then from continnity of the fnnction qoir) we 
will get that this fnnction has a root in the segment [0, 0(h)hs)] C [0, tt]. 

Let us take 


0 (h) Ps) 


TT, for Ir^psl ^ |, 
«, |? 7 P 3 | > i 


If IhPsI ^ 1 then we get |vr?7P3| ^ | and 


go(0(h)P3)) = -cos(7Tr]p3) ^ 0. 


This proves (1) for rj G [—|, |] and the second part of (2). 
If |?7P3| > I then we obtain 0(77,^3) = 2|^F ^ ^ 


?o(0(h)P3)) = -P3sin(0(r7,p3))sin 


-sgn (?7P3 )p3 sin (0 (77, p3)). 


From sin(0(77,p3)) > 0 it follows that 


sgn go(0(h)P3)) = -sgn(?7). 

This proves ( 1 ) for 77 > f. 

Now we need to prove the first part of ( 2 ). The fnnction tq is even, so it is enough to 
prove that ro(p3) ^ tt for ^ ^ P3 ^ 1. 

If P3 = 1 then the equation qo^r) = 0 has the form cos(r(l + 77)) = 0 . The first positive 
root of this equation is 2{i+r]) > tt for — 1 < 77 < —Hence, the statement is satisfied for 
p 3 = 1- 

By contradiction, assnme that there exists pg G [^, 1 ) snch that ro(p3) < tt. From the 
continnity of the fnnction tq it follows that there exists ps G (^3, 1 ) such that Tq^Ps) = n. 
This implies 

?o(P3)7r) = -cos(7r|77|p3) = 0 , 

thus p3 = ^ + 1^, where k E Tj. From — 1<77<—^it follows that for any A; G Z it 
satisfies ps ^ (^, 1 ), hence p^ ^ (pg, 1 ). This contradiction concludes the proof. □ 
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We get the following description of the hrst Maxwell sets (in the image and pre-image 
of the exponential map) that correspond to the symmetries from the group S. 

Corollary 1. ( 1 ) If rj ^ then 

M{S) = {(p,f) e C X R+ I f 

Exp M{S) = P := G SO3 I -y e 7 ^ 0 } ~ RP^ 

is a projective plane that contains axial symmetries. 

{ 2 ) If 7 ] < —\ then 

M(S) = {(p,t) e c X R+ I |p3| < ^, t 

u {(p,t) e C X R+ I \ps\ ^ P3 ^ ±1, t = ^}, 

Exp M{S) = PUL^, 

where 

Lri := {R±es,^ e SO3 I Lf e ( 27 r(l rj), tt]} 

is an interval. 


Proof. The statement on Ai{S) follows directly from Propositions 0 ] and [ 5 l 

To hnd Exp A 4 {S) notice that for any element of SO 3 (i.e., a rotation around some 
axis) a corresponding quaternion has a real part go equal to the cosine of the half of the 
rotation angle. Therefore, if tq is a zero of the function go then the corresponding element 
of SO 3 is an axial symmetry. So we get the component P of Exp Ai{S). 

Let us show that if 77 < — p^ G (— 1 , ~^] U [^, 1 ), r = vr then in the image of the 
exponential map we have Namely, from the formulas for gi,g2 it follows that their 
values are zero for r = tt. Hence, the corresponding transformation is a rotation around 
the axis 63. Next 

<?o(7r) = -cos(7rr7P3). 

This implies that the rotation angle is of the form ± 27 r? 7 P 3 . When ps is in the above 
intervals the rotation angle lies in { 2 'K 7 ], — tt] U [tt, — 2717 ]). If we assume rotations in both 
directions then the set of rotation angles (27r(l + 7 ]), tt] is enough. It is an interval in SO 3 , 
because we need to identify points vr in two half intervals (two rotations by the angle tt 
in different directions are the same transformation). □ 

Let us denote the hrst Maxwell time by 


^max(p) 


2min(7r,ro(p3))/i 

\P\ 


( 11 ) 
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6 Conjugate time 

We recall here the formula for the conjugate time from Lemma 5 [2] and state some 
properties of the conjugate time (Lemma 6 of the same paper). Let tconj(p) = 
be the hrst conjugate time, i.e., the exponential map is degenerate at the point (p, tconj(p)), 
and non-degenerate at any t G (0,fconj(p))- 

Theorem 1 (L. Bates, F. Fasso 0 )- The conjugate time has the following properties: 

(1) // -1 < p ^ 0 then TconiiPs) = vr for any ps e [-1,1]. 

(2) If p > 0 then TconjiPs) is the smallest positive root of the eguation 

1 _ u2 

tanr = -p- -^r, (12) 

1 + wi 

and the ineguality | < rconj(P 3 ) ^ tt zs satisfied. There is the equality only for p^ = ±1. 

(3) The function Tconj : [“Ij 1] —)■ R is smooth. 

(4) The function Xconj increases on the segment [0,1]. 

Proposition 6. The first Maxwell time is less or equal than the conjugate time. 

Figure [3] shows the plots of the hrst Maxwell time and the conjugate time. In the case 
?7 > 0 the plot of the conjugate time belongs to the dashed region. 


Figure 3: The hrst conjugate time and the Maxwell time. 



Proof. From Proposition [5] it follows that if p ^ 0 then the hrst Maxwell time TmaxiPs) ^ tt 
for any p^ G [—1,1]. But from Theorem [1], we get Tconj(Ps) = t. 

From Theorem [H we get that if p > 0 then | < Tconj (ps) ^ t. Let us prove that if 
p > 0 then for any p^ G [—1,1] the inequality To(p 3 ) ^ f is satished. This will complete 
the proof, because from Corollary [T] we get that the hrst Maxwell time is dehned by the 
value of to(p 3 ). 
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The function tq is even, so it is enough to prove the statement when G [0,1]. 

If Ps = 0 then the hrst positive zero of go is equal to | and the statement is satisfied. 
If P 3 = 1 then the first positive zero of go('r) = cos(r(l + p)) is equal to ^ ^ 

when p > 0, and the statement is satisfied. 

By contradiction, assume that there exists pg G (0,1) such that To{p'^) > |. Then 
(since the function tq is continuous) there exists pg G (pg, 1) such that 'ro(pg) = |. This 
means that 

Qo {p3, 1) = -Pssin = 0. 

It follows that Pg = ^, where fc G M and A; < |. 

From formula (j7]) we get the partial derivative of qo{p 3 , t) with respect to the variable 
T at points of the form Xk = (^, f), where fc G Z: 

For /c = 0 we get xo = (0, |) and ^(a^o) < 0. For k = 1 we obtain ^{xi) > 0. Let us 
consider an arc of the graph of the function tq, with ends at the points Xq and Xi (this arc 
exists because of continuity of the function tq). The function t) is smooth, so there 
exists a point (pg, r) on this arc such that ^ is equal to zero at that point. It follows 
that equations ([9]) are satisfied, but this is not possible for pg G [—1,1] as was proved in 
Proposition [3l 

Hence, points of the form pg = ^, where fc G N, do not lie on the segment [—1,1]. 
This contradiction completes the proof. □ 

7 Cut locus 

Let us denote the open set bounded by the closure of the first Maxwell set (corresponding 
to the group S) in the pre-image of the exponential map by 

t/ = {(p, t) G (F X R+ I 0 < t < tmax(p)}- 

Proposition 7. The map Exp : U —?• SOg \ (Exp At (S') U {id}) is a diffeomorphism. 

Proof. We will use the Hadamard global diffeomorphism theorem [10] (a smooth non 
degenerate proper map of smooth arcwise connected and simply connected manifolds of 
the same dimension is a diffeomorphism). 

The sets U and SOg\ (Exp At(S')U{id}) are both three dimensional, arcwise connected 
and simply connected (homeomorphic to a punctured ball). 

From Proposition [6] it follows that the exponential map is non degenerate on U. 
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Let us prove that the map Exp : U —)■ Exp U is proper, i.e., pre-image of a compact 
set K d Exp f/ is a compact set. The set U is bounded, this implies that the pre¬ 
image of K is bounded too. So, it is enough to prove that this pre-image is a closed set. 
By contradiction, assume that there exists a sequence {pn,tn) € Exp“^ K converging to 
{p,t) e I 7 \Exp-i K. 

Then (since the map Exp is continuous) the sequence Exp (p„, G K converges to 
Exp {p,t). 

If {p, t) eU then Exp {p, t) E K (because iL is a compact set). Hence {p, t) G Exp“^ K, 
a contradiction. 

If {p,t) ^ U (i.e., {p,t) lies on the boundary of U) then (since C and A 4 (S) are closed 
sets) t is the Maxwell time or = ±l,t = or t = 0 . Thus, Exp {p,t) lies in 

Exp A 4 (S) or is equal to id. But K is a compact set, so we get a contradiction. □ 

Definition 6. Recall that the cut time tcnt{p) is the time such that a geodesic Exp(p,f) 
is a minimizer for t G [ 0 ,fcut(p)] but is not a minimizer for t > tcut(p)- The point 
Exp(p, tcut(p)) is called a cut point. The cut locus is the set of the cut points. 

Theorem 2. (1) If p ^ | then the cut time is 

^ ^ bsl < 1, 

\P3\ < 2R- 

Theorem 3. (1) If p ^ | then the cut locus is the projective plane 

p = {Rv,tt e SO3 I n G M^, n 7^ 0 } 
consisting of all axial symmetries. 

( 2 ) If p < —I then the cut locus has two components P U Lrj, where 

Lr, = {Rei,±^ e SO3 I p G [ 27 r(l + p), 7 t]} 

is the segment consisting of some rotations around the axis 63 (This axis corresponds to 
the eigenvalue of the metric which is not egual to two others). 

The proofs of Theorems^ and [ 3 ] directly follow from Proposition [71 □ 

The cut locus is the surface of revolution of the figures represented in Fig. 0 ] (in the 
model of SO3 as a ball with antipodal identification of boundary points). 

The wavefronts and the cut locus are shown in Fig. [ 5 ] for p < —I. 


( 2 ) Ifp< 


■I then the cut time is 


2'kI\ 

IpI ’ , 

2/l-ro(p3) 

bl ’ 


for 

for 
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Figure 4 : The cut locus. 








Remark 2. If ^ | then the cut locus consists of the hrst Maxwell points. If r; < — | 

then the cut locus is the union of the hrst Maxwell points and the two conjugate points 

■^e3,±27r(l+7;) • 

Corollary 2. In the Euler case the cut time is equal to vri/Jj", and the cut locus is P. 

Proof. If p = 0 then the cut time is • It is easy to see that if p = 0 then tq = 

and IpI = y/h. 

8 Diameter of SO3 in the Lagrange case 

Now we compute the diameter of a left invariant Riemannian metric on SO3 and describe 
the set of all most distant points from the identity. 
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Theorem 4. (1) The diameter of SO3 is equal to 


27 r/r^^l + ^, for 

TTy/h, for r/e[-i,0], 

for 776(0,+CX)). 

( 2 ) The set of all most distant points from the identity is 

{R±ez,^], for r/G(-l,0), 

P, for 77 = 0 , 

{Pe.TT I e 6 span{ei, 62}}, for 77 G (0,+cx)). 


Proof. The diameter is equal to the maximal value of the cut time. The cut time is 
2Jir^t(P3) ^ rcut(P3) = min(ro(p3),7r). 

Let us compute \p\ as a function of the variable p^. We have 

pI + pI pI . 

h h 

It follows that pi + pi = h — Hence 

H 2 2 I 2 I 2 T 2 

= Pi+P2 + Pz = h- PsV- 

By the substitution p^ = pslp] and expressing |p| we get 

\P\ = 

1. Let 77 G (—1, — |). The hrst Maxwell time that corresponds to the rotations around 
the vertical axis is equal to This time (as a function of the variable ps) increases on 
the segment [—1, 0] and decreases on the segment [0,1]. Moreover, this time is greater or 
equal then the cut time on the segment and outside this segment they are equal. 

This implies that the cut time has the maximum value at ps = This maximum 

value is equal to 1 + ^. The corresponding points with maximal distance to the 

identity are {P±e 3 , 7 r}- 

2 . Let 77 ^ In this case the cut locus is P. It is clear that the set of all most 
distant points from the identity is a subset of the cut locus. 

Let us consider the set of curves of the form {Pp,t | t G [0,7r]}. These curves connect 
the identity with a point of the cut locus. Let us compute the length of such curve: 

hPi + hPl + hPldt = Ti\/h 
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It is obvious that this length as a function of the variable is monotonic on the segments 
[— 1 , 0 ] and [ 0 , 1 ]. Furthermore, this length is greater or equal than the cut time. And the 
equality is satisfied for p^ = ±1 or ps = 0 . It follows that the diameter of SO3 is equal to 
the maximum value of the cut time at these points 

rcut(O) = rcut(±l) = 2(1 + ^)- 
The corresponding cut times are 

tcut(O) '\/11 ) tcut(Tl) ^ • 

2 a. If p G [—|, 0 ] then the maximum value of the cut time is ti'/Ts. If p 7^ 0 then the 
set of all most distant points is {i?±e3,7r}, and if p = 0 then the set of all most distant 
points is the projective plane P. 

2 b. If p G ( 0 , +00) then the maximum value of the cut time is The set of all 

most distant points is a circle {Ry^n I |n| = 1 , u G span{ei, 62}}. □ 

9 Left invariant Riemannian problem on SU 2 
in the Lagrange case 

Let us consider the left invariant Riemannian problem on SU2 in the case of two equal 
eigenvalues of a metric. In the case of p > 0 the result of T. Sakai [ 3 ] is that the cut locus 
is a two dimensional disk. In the case of p < 0 there is a conjecture (M. Berger | 3 ]) that 
the cut locus is a segment. If p = 0 then this segment becomes the point id. We will give 
the proof of this conjecture. 

Let us use the same method, as in case of SO3, to hnd the cut locus. First, notice that 
the exponential map can be written by the same formulas (jl]). Secondly, the symmetry 
group of the exponential map is the same as in the case of SO3. The difference is that 
the set Exp Aio is not a Maxwell set on SU2. In the SO3 case there are two geodesics 
that come to a point of this set at the same time, but when we lift them to SU2 these 
geodesics are on different leaves of the covering at that time. 

But the conjugate locus and the conjugate time have the same description as in the 
SO3 case. Therefore, for application of the Hadamard global diffeomorphism theorem, 
the main question is the comparison of the Maxwell time for the Maxwell strata Exp Ad 12 
and Exp and the conjugate time. The answer is in Propositions 151 and [Till A technical 
Proposition [ 9 ] is needed to prove Proposition [101 

Proposition 8. If —1 < rj ^ 0 then for any p^ G [—1,1] \ {0} the inequality T^{pz) A tt 
is satisfied. 
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Proof. It is enough to prove the statement for pa G (0,1], because the function ra is even. 

It is clear that ra(l) = and the statement is satished for pa = 1- By contradiction, 
assume that there exists pg G (0,1) such that < tt. It follows that there is pa ^ 

(p'a, 1) such that ra(pa) = tt (because the function ra is continuous). By substituting this 
value to the formula of ga, we get 


— sin(7r?7P3) = 0. 


Hence, ps = j^, G Z \ {0}, it follows that pa ^ (0,1). We get a contradiction. □ 

Proposition 9. // p > 0 then for any £ > 0 there exists 0 < |pa| < £ such that the 
inequality ra(pa) ^ 'rconj(O) is satisfied. 

Proof. From dehnition of the conjugate time we see that the point (0,rconj(0)) is the 
critical point of the function q 3 {p 3 ,T). 

It is easy to compute the second partial derivatives of this function: 


d\3 

dpi 


—TT] sin(rppa) (rp cos r + 2 sin r) — r^p^pa sin r cos(rppa), 


d\3 

dTdp3 


((l + p)cosr—(1 + ppa) r p sin r) cos ( Trip 3 ) — ppa sin (rppa )((H-p)rcosr + 2sinr), 


dr"^ 


(1 + 2 pp 3 + p Pa) COSTsin(rpp 3 ) — p 3 (l + 2p + p Pa) sinr cos(rpp 3 ). 


Hence, the Hesse matrix of the function ga at the point (0, rconj(O)) is (since sinrconj(O) = 
-'rconj(0)pcosrconj(0) by Theorem [I]) 


( 0 (1 + p + (O)h^) cos Tconj (0) 

V (1 + h + Tonj (O)h^) cos Tconj (0) 0 

Because of | < Tconj (0) < tt (Theorem [T]), the nondiagonal elements of this matrix are 
nonzero. So, the point (0, Tconj(0)) is a saddle point of the function ga. Thus, two isolines 
of the function ga intersect transversally at the point (0, Tconj(O)). One is the line pa = 0 
and the second consists of points (pa, t) where t is some positive root of the equation 
g 3 (p 3 ,T) = 0. If it is not the smallest positive root then there exists e > 0 such that for 
all 0 < Ipal < e we have T 3 (p 3 ) < Tconj(0). 

Consider now the case when the second isoline consists of points (p 3 ,T 3 (p 3 )). By 
contradiction, assume that there exists a punctured neighborhood of zero where T 3 (p 3 ) > 
'Tconj (0). Then in this punctured neighborhood we have the inequality 



T'3(P3)h + tanT3(p3) > 0. 


(13) 
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Let us calculate the derivative of the function 


dTs 

dp3 


(Ps) 


dqs j dq^ Trj cos r cos^rrjp^) + sin r cos^rrip^) — rpp^ sin r sin^rpps) 

dps/ dr —{1 + pPD sin Tsm{Tpp 3 ) + ps{l + p) COST cos{Tpp 3 ) 


Note that due to continuity of the function ts we can choose a punctured neighborhood 
of zero such that there holds the inequality cos r cos(r?7P3) 7^ 0 . After dividing the terms 
of the fraction by cos r cos^rpps) we get 


dTs 

dps 


(Pa) 


dqs / dqs rp + tan r — rpps tan r inni^pp^) 

dps/ dr -(1 + 77^3) tanrtan(r?7P3)+p3(l + r/)’ 


If r = T^{p^) then 


tan(r77P3) + ps tan r = 0. 


This implies 


dr^ rp + tan r + rpp^ tan^ r 

dps (1 + VPs)p3 tan^ r + ^ 3(1 + 77) ‘ 


( 14 ) 


From inequality (IT^ and 77 > 0 it follows that the sign of this derivative is opposite to 
the sign of 773. This means that in some punctured neighborhood of zero the function 
increases for 773 < 0 and decreases for 773 > 0 , also Jim T^i^ps) = Tconi(O), and by assumption 

P3^0 

'I'siPs) > T'conj(O) in this neighborhood. The contradiction completes the proof. □ 


Proposition 10. If p > 0 then for any p^ G [—1,1] \ {0} the inequality r^^ps) ^ Tconj(P 3 ) 
is satisfied. 


Proof. Let us assume that 773 G ( 0 , 1 ]. Because the functions T 3 and Tconj are even, this is 
enough to prove the proposition. 

If P3 = 1 then we get r3(l) = < vr = Tconj( 1 ), he., the statement of the proposition 

is satished. 

By contradiction, assume that there is pg G ( 0 , 1 ) such that rfip'fi) > Tconj(Ps)- From 
Proposition [3 it follows that there exists 773 G (0,^3) such that rfipfi) < Tconj(P 3 )- 

This implies that on the segment [773, 1 ] there exist at least two points with the equal 
values of the functions Tg and Tconj (because the functions T3 and Tconj are continuous). 

Notice that if ps G ( 0 , 1 ] is such that rfipfij = Tconj (7^3) then ^{ps) < 0 . This follows 
from formula flTT)) . because Tconj( 0 ) < Tfip^) < n (the function Tconj is increasing), hence 

^3(P3)P + tanT3(p3) > 0. 

So, at all these points (there are at least two such points) the derivative of the function 
T3 is negative. We get a contradiction. □ 


Theorem 5. The cut locus on SU 2 is 
(1) the segment 


Tr, := {-cosiupps) - sm{Tipps)k G H | p 3 G [-1,1]} 
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for —1 < ?7 ^ 0 (if rj = 0 then becomes the point { — 
(2) (T. Sakai J^) the disk 


Dn := {qo{p,n{p3)) + qiip, T^{pz))i + q2{p,T3{p3))j G M | p G C}, 
bounded by the circle of conjugate points 


S := {cosrconj(O) + sinrconj(0)(icos(^ + j sin^j) G H | G [0,27r]}, 


for T] > 0. 

Proof. From Propositions [ 8 ] and [10] and from Theorem [1] it follows that the exponential 
map is non degenerate on the open set 


{{p,t) I p G C, 0 < f < 


for p ^ 0 and on the open set 



for ?7 > 0. The image and the pre-image of the exponential map are arcwise connected and 
simply connected. The proof of properness of the exponential map is qnite the same as in 
the SO 3 case. So by the Hadamard theorem HDl the exponential map is a diffeomorphism 
on these open sets. 

Tims the cnt locns is the image by the exponential map of the following sets, respec¬ 
tively: 




It is easy to see that images of these sets are and respectively. 


□ 


10 Connection with the sub-Riemannian problem 

on SO 3 


By identifying the Lie algebra 503 with the space of imaginary qnaternions we consider 
the decomposition 


503 = !©p, 


(15) 


where t = Mfc, and p = Mi © Mj. 
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Let A be the distribution on the Lie group SO 3 that is generated by left shifts of the 
subspace p. Let us endow A with a positive dehnite quadratic form rg(v) = (g~^v, g~^v), 
where g G SO 3 , v E Ag = gp, and (•, •) is the Killing form. Let Xi,X 2 be vector helds 
that he in the distribution A and form at any point g G SO 3 an orthonormal frame of Ag 
with respect to the form Vg. 

Let us consider the following left invariant sub-Riemannian problem: 



(Actually there are many other left invariant sub-Riemannian structures dehned by the 
distribution and a Riemannian structure on it. Their classihcation for three dimensional 
Lie groups was obtained by A. A. Agrachev and D. Barilari cao 

Theorem 6. For the left invariant Riemannian problem on SO3 in the Lagrange case the 
following objects converge to the corresponding objects of the left invariant sub-Riemannian 
problem on SO3, defined by decomposition 03 and the Killing form, as J3 —)• oo.- 

(1) the parametrization of the sub-Riemannian geodesics, 

(2) the conjugate time, 

( 3 ) the conjugate locus, 

( 4 ) the cut time, 

( 5 ) the cut locus. 

The cut loci of the sub-Riemannian and Riemannian metrics (in the case rj < —i) are 
shown in Fig. [6] (the surfaces of revolution of the represented hgures). 

Proof. ( 1 ) Parametrization of Riemannian geodesics on SO3 has the form 



where p G SO3 and p = pR -\- P2J + Psk is the corresponding imaginary quaternion. If 
1] ^ —1 (this is equivalent to I3 ^ 00) then we get 



This coincides with a well-known parametrization of sub-Riemannian geodesics (see proof 
in V. Jurdjevic’s book [ 13 ]) 


g{t) = exp{t{Ap + Ak)) exp{-tAk), 
where Ak El, Ap Ep, rid(Ap) = 1. 
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Figure 6: The cut loci in the sub-Riemannian and Riemannian cases. 




In paper [7] U. Boscain and F. Rossi consider an initial co-vector 

- cos 6i -\ — sin 0? -|- ck. 

2 2 

In our work the initial co-vector p has a form 

. , , pI , pI , pI 1 

Pl^+P23+P3k, — + — + — = 1. 

^2 ^3 

It is easy to see that in the Lagrange case (Ji = I 2 ) 

bP = h- vpI- 

Hence, \p\ —)■ s/l + as r/ —)■ — 1 and R = 1. 

(2) In the paper by L. Bates and F. Fasso [2] (Lemma 5) there is an equation for 
the hrst conjugate time for the Riemannian problem on SO3 in the Lagrange case. If 
— 1 < 77 R 0 then the hrst conjngate time is equal to tt. If 77 > 0 then the hrst conjugate 
time is the smallest positive root of the equation 


tanr = —p 


1 - pi 


T. 


1 + ppI 

As 77 —)■ —1 we get the equation for the conjugate time in the sub-Riemannian problem 



















where r = ^|p| —)■ \/l + | as —)■ —1. The hrst conjugate time corresponds to the 

root of the hrst factor of this equation [ 2 ] . 

(3) If -1 < r] < 0 then in the Riemannian case the conjugate locus is a segment of 
the length 47r\ri\ or the circle (see Proposition 2 in |2]) 

Srj = {exp(±(p/c) I (p e [-27r|r7|,27r|r7|]}. 

If —1 < T] ^ i then it is a circle, i.e., as —)■ — 1 it converges to the conjugate locus in 
the sub-Riemannian case [7]. 

(5) If ?7 —)■ —1 then the cut locus component of the Riemannian problem 

Lr^ = {exp(±(p/c) I if e [27r(l -h //), vr]} 

converges to the circle = {exp{ipk) | G M} (the cut locus component of the sub- 
Riemannian problem). The rest part P (the set of all axial symmetries) of the cut locus 
is the same for both problems. 


(4) Let us calculate the time at which geodesics reach every point of the cut locus in 
the Riemannian case. 

1. Consider the component of the cut locus Let a = qo + q^k G El define a point of 
The cut time of Riemannian problem for this component is defined by = tt. We 
get 

f - cos 7Tp3 = qo, 

I sin7rp3 = g3. 

It follows that arga = tt — vrps. Hence 




h 

i + pIp 


__ 

7r2(l + p) — 2'Kpaxga + {aigaYp 


Thus the cut time for this component is 


2'jiJ 

tent = -r-r = 2A/iT\/vr2(l + p) - 27rargr7a {aigayp. 

\P\ 


As p ^ —1 and R = 1 the cut time converges to 2-^arga(27r — arga), and this is the cut 
time for the sub-Riemannian problem. (This easily follows from Theorem 2 of [7], which 
contains formulas for the sub-Riemannian distance.) 


2. Consider now the component of the cut locus P. Let us calculate the time at which 
the geodesics reach the point Il{qii + q 2 j + q^k) G P. Denote a = q^i, (3 = qi + q 2 i & H. 
From parametrization of Riemannian geodesics we have 

siuT^pl+pl = \(3\. 
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From 


pI+pI = 


\/pI+pI \/h- {p + l)p\ 


\p\ 


it follows the equation for pa 

sinr-^/i - (p + l)p: 


\P\ 


\P\ 


= \j3\ = 


The equation for the cut time is 

cos r cos(r?7P3) — p3sinrsin(r?7P3) = 0, 

it is equivalent to 

pgtanr = cot(rpp3) = tan Q - rppgj , 


or 


TT 


rpps + arctan(p3 tanr) = — = arga. 


As ?7 —>■ — 1 we get 


IpI \/l + c2, p 3 


r —)■ 


\/l + 


■ 2 

hence we have the formulas of the sub-Riemannian distance [ 7 ] 


ct , c , \/l + . 

+ arctan(^q-^tan(-^-)) = arga, 


sm 


• Vl+c^f' 


\/l + 




a 


□ 


Conclusion 

The general left invariant Riemannian problem (with Ii < I2 < I3) is much more compli¬ 
cated that the Lagrange case considered in this paper. In the general case the problem 
has no rotational symmetry, and geodesics are parameterized by non-elementary functions 
(elliptic integrals). Although, we believe that some results on optimality of geodesics can 
be obtained by the approach used in this work. 
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